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The optical conductivity (j{uj) of a doped two-leg t-J ladder is calculated for an electric field 
polarized parallel to the legs of the ladder. The conductivity has a Drude weight proportional to 
the hole doping and an apparent threshold for absorption (a pseudo gap) which may be associated 
with the energy to break a pair. This pseudogap in a (a;) is present even though the pairs have a 
modified d^2_y2-like wave function because the geometry of the ladder leads to quasi-particle states 
which probe the gap along an antinode. 
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Theoretical predictions^ of the existence of a spin 
gap in a two-leg S — 1/2 antiferromagnetic ladder 
have now been observed in spin Sjjsceptibility measure- 
ments of vanadjd pyrophosphates {VO)2P20^ and the 
cuprate materialcl SrCu203. These materials are insu- 
lating so that the charge excitations are gapped as well 
as the spin excitations. Recently, Hiroi and TakanocI 
have succeeded in doping a two-leg ladder compound 
Lai^xSr^u02.5, andrnQbserved an insjtikfctor to metal 
transitionla. NumericalHlJ'El and analyticBll3 calculations 
have found that hole-doped two-leg t-J ladders have a 
gap to spin excitations but no gap to charge excitations. 
These calculations show that the doped holes form mod- 
ified 'c?a;2_j,2-like' pairs and that the pair-field correla- 
tions have a power-law decay. In addition the pairs are 
expected to exhibit power-law 'ikf' charge-density cor- 
relations. 

The absence of a charge gap implies the existence of a 
finite Drude weight which is expected to be proportional 
to the hole doping at small doping levels. In addition, the 
possibility of breaking up the pairs suggests that in spite 
of the vanishing of the charge gap reflected in the col- 
lective motion of the pairs, there should be a pseudogap 
in the frequency dependence of the charge excitations 
associated with breaking the pairs (we note that Gia- 
marchi et qlj have predicted oj^ behaviour for the low 
frequency conductivity of a Luttinger Liquid, although 
direct calculations of this behaviour are as yet unob- 
served). Here we seek a further understanding of the 
charge response of the doped two-leg t-J ladders by ex- 
amining the frequency-dependent conductivity cr(Lu) for 
an E-field polarized along the legs of the ladder. Lanczos 
calculations for doped 2x5 and 2 x 10 t-J ladder have 
been carried out and the behaviour of a{Lu) has been in- 
vestigated. 

The t-J hamiltonian on a 2 x L ladder is defined as 
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where f5 (=1,2) labels the two legs of the ladder (oriented 
along the s-axis) and j is a rung index (j=l,...,L). Sj^p 



and 



are electron spin and creation operators and 



the Gutzwiller projector Pq excludes configurations with 
doubly occupied sites. In the materials of interest, the 
exchange coupling J along the legs is nearly the same 
as the exchange coupling J' across a rung, and similarly 
the hopping t along the legs is close to the rung hopping 
strength t'] therefore, in the following calculations we will 
work with the isotropic system, J = J' and t = t' . 

For an electric field polarized along the a;— axis (paral- 
lel to the legs) , the optical conductivity can be expressed 
as 



a{uj) = 2ttD5{uj) 
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where the final term is usually defined as the regular part 
of the conductivity {a^'^^{uj)). D is the charge stiffness 
and jx is the paramagnetic current operator 
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where we choose the charge e to equal unity. = 2i is 
the number of sites and M is the number of holes in the 
ladder (assumed to be even throughout). The states \(t>^) 
in Eq. ^ are the energy eigenstates of the hamiltonian 
7i with corresponding energies 

En'^ is the ground 

state with energy E^^ . 

The Drude weight ctq = 27r£' may be evaluated in two 
ways. Firstly, as explained in a previous publicationB, 
the charge stiffness D can be calculated by considering 
the curvature of the ground state energy level [E^^] as 
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a function of the flux <I> (in units of the flux quantum 
<I>o = hc/e) threaded through the two chain system 



D = 
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Secondly, (Tq may be calculated using the sum rulfJ^^'t^' 
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where the Kinetic Energy operator Txx along the chain 
direction consists of the term proportional to t in the 
hamiltonian (|l|). Substituting Eq. ^into Eq. we have 
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where the second term is equal to {N/n) a^'^^ {uj)dLO. 

In Fig. 1^ we plot N x f7'"'^^(a;) for the ladder system 
with various values of the ratio J/t: The results were ob- 
tained using the Lanczos approach and data from both 
the 2x5 and 2 x 10 system are shown. In order to ob- 
tain the absolute ground state of the system, we have 
chosen boundary conditions to form a closed shell in the 
non-interacting Fermi sea (obtained by turning off the 
interactions J and J'): specifically this corresponds to 
anti-periodic boundary conditions for n < 0.5 and pe- 
riodic boundary conditions for n > 0.5. We have also 
considered the parity of the states under a reflection in 
the symmetry axis of the ladder along the direction of 
the chains (even {R^ = 1) or odd {Rx = — !))■ 

In Fig. 0a we plot the Drude weight ctq as a function 
of electron density for various values of the ratio J/t; we 
note that both of the methods above were used (Eqs. ^ 
and H) and the results were in excellent agreement. In 
Fig. we plot the integrated finite frequency conduc- 
tivity 2 dujcr^'^^ {cj) as a function of electron density 
for various values of the ratio J/t; this quantity can be 
directly related to ctq and to tt times the kinetic energy 
per site through the sum rule (Eqs. || and||). 

There are several important features of Figs. |l| and ^ 
we should mention. We have checked that as the system 
size is increased, the ratio of weight at finite frequency to 
the total conductivity remains effectively constant; the 
difference in the number of delta peaks as the system 
size is changed is merely a finite size effect. We have also 
checked that the moments of the finite frequency distri- 
bution are effectively unchanged with increasing system 
size. A 'pseudogap' appears in the optical conductiv- 
ity, below which there appears to be no weight; we dis- 
cuss this feature in some detail below; note that as the 
ratio J/t is increased, the 'pseudogap' increases. Sec- 
ondly, as can be seen by comparing Figs. ||a and b, 
the weight of the conductivity at finite frequency is ex- 
tremely small as compared to that at zero frequency (the 
Drude weight). This behaviour is reminiscent of the one- 
dimensional t-J model and in contrast to the larger finite 



frequency weight found in the two-dimensional systemS. 
The Drude weight is proportional to the hole-doping at 
low doping levels, and at low densities is effectively inde- 
pendent of J as we would expect for a spin-charge sepa- 
rated state. Note however, that the small finite frequency 
absorption has a larger magnitude together with a larger 
and opposite J dependence for n > 0.5. These differences 
might be related to the nature of the ground state, a 
Tomonaga Luttinger liquid at low electron density and a 
Luther-Emery like liquid with a spin gap at small hole 
density. 

In picturing the intermediate states which enter Eq. |^ 
for the conductivity, it is convenient to consider the limit 
in which the rung exchange J' is large compared to both 
J and the hopping t{=t'). In this case, when two holes 
are doped into the antiferromagnetic ladder they will go 
onto the same rung in order to minimize the number of 
broken singlet rung bonds. The binding energy of this 
state is approximately J' —2t—2t' , reflecting J' the energy 
required to break another rung singlet minus the kinetic 
energy gain for two separate holes. In the limit of large 
J', the doped system can therefore be thought of in terms 
of hole pairs moving on a lattice of rung singlets with an 
effective pair transfer matrix element of order —2t^/J'. 
The fact that only one pair of holes can occupy a given 
rung leads to a low-energy description in terms of hard- 
core charge 2e bosons. 

The coherent propagation of pairs leads to a finite 
Drude weight. In addition, when w exceeds an energy 
of the order of the pair binding energy Eb, it is possible 
for the system to absorb a photon and excite a state con- 
taining two quasi-particles (each with charge e and spin 
1/2). The two quasi-particle state which enters Eq. |^ for 
a{u!) has S=0 and a center of mass momentum equal to 
zero. There is in fact a continuum of two quasi-particle 
scattering states above the threshold energy for breaking 
a pair in the infinite ladder. As discussed by Tsunetsugu 
et aB, a similar two quasi-particle state with S = 1 de- 
termines the spin gap in the doped two-leg ladder. An 
estimate of the pair binding energy can be obtained by 
considering the ground state energy of a system with an 
even number of holes M and the corresponding energy 
if a hole with a momemtum k along the chains is either 
added or removed. Such excitations are infact the quanti- 
ties we would expect to measure in the spectral functions, 

and A{k, c) = E„ I {<l>^f^'\ck,s m?S{uj + E^+^ - E^) 
where the first expression relates to the addition of an 
electron and the second relates to the addition of a hole. 
Indeed, sharp features in the spectral function should ap- 
pear at energies e^^^ — ±(i?o^ - £'g^*^(fc, i?a;)) where 
we can obtain two different values depending whether 
look at bonding or antibonding excitations {Rx = ±1)11 

In order to consider in more detail the quasi-particle 
excitations we define the (average) chemical potential as 
/i=l/2(£'(J^+^ - E^'^^) where E^ relates now to the ab- 
solute ground state energy for P holes irrespective of the 
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quantum numbers k and Rx- In Fig. ^ we show the 
dispersions — /i for the 2 x 10 ladder at n = 0.8 

and J/t ~ 1.0 with positive energies corresponding to 
the case of adding an electron to the 4 hole ground state 
and negative energies corresponding to adding a hole to 
the 4 hole case. A finite gap separating the electron and 
hole quasi particle spectra can be clearly seen on the fig- 
ure both for the bonding or antibonding sectors. Indeed, 
these gaps are related to the binding energies of a pair, 

Eb{R,) = E^'-\Rx) + E'^'+\R,) - 2El', (7) 

i.e. the energy required to create an electron-hole exci- 
tation with Rx = ±1. These two processes correspond in 
Fig. H to transitions involving the smallest energy exci- 
tations as indicated by the arrows. 

Since these transitions conserve momentum (not quite 
exactly for R^ = —1) they can be considered as 'verti- 
cal' optical transitions and should be seen in the opti- 
cal absorption (H). In Fig. ^ we compare the binding 
energies, the pseudogap obtained from the optical con- 
ductivity data and the two spin gaps (corresponding to 
bonding and anti-bonding spin excitations) at n = 0.8 
(i.e Nh = 4) for various ratios of the parameter J/t in 
the 2 X 10 ladder. As expected the threshold for absorp- 
tion is of the order of the pair binding energy and also 
provides a measure of the spin gap in the R^ — 1 spin 
excitation spectrum. This adds credit to our picture of 
the pseudogap being a threshold for unbinding pairs. It 
is interesting to note that other lower energy gapped spin 
excitations exist (with R^ = — 1) which probably do not 
involve pair breaking. 

In conclusion, we find that cr(w) for the two leg ladder 
has a Drude contribution consistent with a coherent pair 
motion and an absorption onset which we associate with 
the creation of two quasi-particles. It is interesting to 
note that for the ladder there is a pseudo-gap in (j{u>) 
even though the pairs have a ci3,2_y2-like wavefunction. 
Basically, one is confined by the geometry to probe the 
gap along the antinodc so that cr(w) has a finite threshold. 
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A calculation of the relative weights of the quasi-particle 
excitations is left as a future study. 

FIG. 1. TVx a''^%uj) for 2 X 5 and 2 x 10 ladders at electron 
density n = 0.8. The different plots correspond to a) J/t=0.5 
b)J/t=1.0 c)J/t=1.5 

FIG. 2. a) The Drude weight (tq, and b) The sum of the fi- 
nite frequncy conductivity 2 <^{^), versus electron density 
for various values of J/t. These two quantities are related by 
the sum rule (see text). 

FIG. 3. Energy dispersion for the 2 x 10 ladder system 
{J/t = 1.0) with 3 holes (upper curves) and 5 holes (lower 
curves); the results for both = 1 and = — 1 are shown. 
Momentum is in units of tt. The arrows represent the en- 
ergy required to create an electron-hole excitation with either 
Rx = ±1 (see text). 

FIG. 4. A comparison of the the binding energy, the pseu- 
dogap (obtained from the optical conductivity data) and the 
spin gaps. These quantities are defined in the text. 
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